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ABSTRACT 

Precise astronomical spectroscopic analyses routinely assume that individual pixels in charge- 
coupled devices (CCDs) have uniform sensitivity to photons. Intra-pixel sensitivity (IPS) vari- 
ations may already cause small systematic errors in, for example, studies of extra-solar planets 
via stellar radial velocities and cosmological variability in fundamental constants via quasar 
spectroscopy, but future experiments requiring velocity precisions approaching ~1 cm s"' will 
be more strongly affected. Laser frequency combs have been shown to provide highly precise 
wavelength calibration for astronomical spectrographs, but here we show that they can also be 
used to measure IPS variations in astronomical CCDs in situ. We successfully tested a laser 
frequency comb system on the Ultra-High Resolution Facility spectrograph at the Anglo- 
Australian Telescope. By modelling the 2-dimensional comb signal recorded in a single CCD 
exposure, we find that the average IPS deviates by < 8 per cent if it is assumed to vary sym- 
metrically about the pixel centre. We also demonstrate that series of comb exposures with 
absolutely known offsets between them can yield tighter constraints on symmetric IPS varia- 
tions from ~100 pixels. We discuss measurement of asymmetric IPS variations and absolute 
wavelength calibration of astronomical spectrographs and CCDs using frequency combs. 

Key words: instrumentation: spectrographs - instrumentation: detectors - methods: labora- 
tory - techniques: spectroscopic 



1 INTRODUCTION 

A growing variety of measurements in astrophysics rely on, or 
are enabled by, precise astronomical spectroscopy. For example, 
searches for extra-solar planets via the radial velocity method (e.g. 
[Mayor & Queloz|1995[rMarcy & Butler|1996^ rely on the repro- 
ducibility of stellar radial velocities over time-scales up to several 
years. Night-to-night reproducibility better than 1ms"' has been 
achieved when integrated over spectral ranges of ~ 1000-3000- A 
(e.g. |Lovis et al.|2006j ). Another example is the use of relative ve- 
locity shifts between ionic and molecular transitions in highly red- 
shifted quasar absorption systems to constrain cosmological varia- 
tions in fundamental constants, such as the fine-structure constant 
(e.g. Ba hcall et al.|1 967; Webb et al. 19991 and proton-to-electron 
mass ratio (e.g. Varshalovich & Levshakov ,,19931 |lvanchik et al.| 
|2005[ l. Velocity precisions averaged over several transitions in hun- 
dreds of absorption systems reach ~20ms"' i Murphy et al.|2063 
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Much more precise and accurate spectroscopy would be re- 
quired if the acceleration of the Universal expansion is to be mea- 
sured by monitoring the slowly changing redshifts of Lyman-or for- 
est absorption lines in quasar spectra (e.g. |Sandage||1962[ |Loeb| 
|1998[ l. This is one possible aim of the CODEX spectrograph pro- 
posed for the European Extremely Large Telescope (e.g. |Pasquini| 
|et al. 2006 , Liske et al. 2008 1. Two epochs of spectra are required, 
taken decades apart, with velocity precision approaching 1 cms"' 
integrated over ~ 1000-3000 A ranges. The calibration of individual 
quasar exposures would therefore need to be reproducible to within 
~1 cms"' over several decades. Many optical and mechanical ele- 
ments of the telescope and spectrograph may need replacing on 
such a time-scale, or observations from other telescopes may need 
to be employed to achieve the required signal-to-noise ratio (SNR). 
Thus, it would be preferable to absolutely calibrate the wavelength 
scale of each quasar exposure. 



Laser frequency combs (LFCs; e.g. jReichert et al.|1999[pones| 
let al.|2000||Udem et al.|2002) are a promising candidate for such 
highly precise, and possibly absolute, wavelength calibration of 
astronomical spectrographs ( [Murphy et al.||2007[ [Schmidt et al.| 
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|2008[ >. The train of pulses that are output from a frequency comb 
appear in frequency space as a spectrum of Unes "modes" uni- 
formly spaced by the pulse repetition rate. This repetition rate can 
readily be synchronised with an absolute radio frequency refer- 
ence, e.g. a global positioning system (GPS) signal or atomic clock. 
Thus, an absolute, very dense and perfectly regular calibration 
spectrum could be used to calibrate astronomical spectrographs. 
[Murphy et al.| \2QQ1) demonstrated that, in principle, recording 
such a frequency comb signal on a conventional charge coupled 
device (CCD) typical of modern astronomical spectrographs would 
provide enough spectral information to reach the ~lcms"' cali- 
bration precision required for measuring real-time drifts in the red- 
shifts of Lyman-Q- absorption lines in quasar spectra. 

Prototype frequency comb calibration systems are now be- 
ing tested at various observatories. The first LFC demonstration on 
an astronomical spectrograph by |Steinmetz et alT] (21308) achieved 
state-of-the-art calibration precision of ~9ms"' in near infra-red 
(IR, 1.5 (i.m) astronomical spectroscopy. |Wilken et al.| ( |2010j l used a 
frequency-doubled Yb-doped fibre laser to calibrate a single echelle 
order at ~500 nm of a highly stable vacuum spectrograph (HARPS) 
with 15 ms"' accuracy, precision and reproducibility over several- 
hour time-scales. [Benedick et al.|(2010| l calibrated the TRES spec- 
trograph with an accuracy better than 1ms"' at ~400nm using a 
tunable Ti:Sapphire LFC. While these tests demonstrate the main 
advantage of LFC calibration - the high density of high-SNR, 
equally-spaced modes - some challenges discussed by ^ Murphy^ 
[eTaTlpOOTl i still remain. For example, the above tests employed 
LFCs covering only small wavelength ranges (typically ~10nm). 
Simultaneously and uniformly illuminating the entire optical wave- 
length range covered by modem echelle spectrographs is highly 
desirable. 

Another important advantage of LFC calibration is that the 
high density of recorded modes and the absolute frequency scale 
should enable direct characterization of many, if not all, instrumen- 
tal systematic errors ( Murphy et al. 2007 ). For example, in calibrat- 
mg a sing le echelle order of HARPS, Wilken etaT] ( |20T0? found 
discontinuities in the calibration curve (i.e. wavelength vs. spectral 
pixel) at regular intervals of 512 spectral pixels. They attributed 
this to imperfections in CCD pixel size, shape and/or sensitivity 
imparted to every 512th pixel during the manufacturing process of 
the CCD mask. Even subtler problems in the CCD or spectrograph 
optics are likely to limit spectroscopic accuracy unless they are first 
detected, measured and then removed with LFCs. 

One such problem is possible intra-pixel sensitivity (IPS) vari- 
ations. Individual CCD pixels can show large spatial variations in 
sensitivity to photons depending on the manufacturing process and 
depth of the depletion layer jToyozumi & Ashley 2005 1. Pixels can 
also vary in size or shape and may contain small defects. Elec- 
tron counts may also be registered in a pixel when a photon enters 
near the edge of a neighbouring one, depending on incidence angle. 
These, and many other effects, lead to apparent IPS variations. The 
development of 'deep depletion' CCDs has mitigated many of these 
effects, but even moderate IPS variations compared to those found 
by [Toyozumi & Ashley] (2005 [ > might affect measurements of, for 
example, possible variations in the fundamental constants. Thus, it 
is vital to accurately measure the detailed behaviour of CCD chips 
- even of individual pixels - if photon-limited spectroscopic preci- 
sion is to be demonstrated. 

This paper reports a new demonstration of an optical LFC us- 
ing a very high resolution astronomical spectrograph with the aim 
of measuring in situ any IPS variations in a modern CCD. Section|2] 
explains our experimental setup. Section|3]sets out the basis for our 



analysis to extract the IPS variations from the comb exposures. Sec- 
tion |4] contains our main analysis and results. Section |5] describes 
how series of LFC exposures might allow absolute wavelength cal- 
ibration of CCDs. We summarize our main results and conclude in 
Section|6] 



2 EXPERIMENTAL SETUP 

To demonstrate that the intra-pixel sensitivity (IPS) variations in 
an astronomical spectrograph's CCD can be measured in situ, we 
tested an LFC at the Anglo-Australian Telescope (AAT), Australia, 
on the 9th April 2010. The Ultra-High-Resolution Facility (UHRF) 
spectrograph has a very high resolution mode, with resolving power 
R = ^/FWHM « 10^ (where FWHM is resolution element's full 
width at half maximum), which means a relatively simple LFC can 
be employed (Section [2.2| l. 

Our experimental setup is illustrated in Fig. [T| After summa- 
rizing the basics of LFCs immediately below, we describe the im- 
portant elements of our particular setup, including our LFC, the in- 
terface with the UHRF spectrograph and the CCD whose IPS varia- 
tions we seek to measure. We also describe the data reduction steps 
taken before the analyses described in subsequent sections. 

2.1 Laser frequency comb basics 

A laser will generate light at discrete frequency modes throughout 
the spectral bandwidth of its gain material as long as the optical 
gain exceeds the losses for a full round-trip in the laser cavity at 
each of the mode frequencies. The frequency of a mode is deter- 
mined by ensuring that there is a whole number of wavelengths 
in one cavity round trip. Conventionally, much effort is expended 
in ensuring that only a single mode can circulate at any time in 
the cavity; however, for mode-locked lasers, one attempts to have 
many modes circulating simultaneously. In addition, by including 
an intensity-dependent loss in the cavity it is possible to bring about 
a mode-locked state in which the laser's output energy is delivered 
in a single, intense and temporally very short pulse. In this condi- 
tion each mode maintains a fixed and well-defined phase relation- 
ship with all other modes. The frequency spacing between all pairs 
of adjacent modes, called the repetition frequency, is exactly the 
same and equal to the reciprocal of the round-trip pulse circulation 
time in the cavity. Further details of the operation of mode-locked 
lasers can be found in Cundiff & Ye|p003t . 

The output spectrum of the mode-locked laser appears as a 
comb of modes where the «th mode has a frequency, given by 

f„ = n/icp + ,/cEO , (1) 

where n is typically of the order of lO"". The repetition frequency, 
/rep, is the mode spacing between the 'teeth' of the comb. The 
carrier-envelope offset frequency (/ceo) represents the offset of 
each mode from an exact harmonic of the repetition frequency, 
while in the time domain it expresses the rate at which the relative 
phase of the carrier and pulse envelope evolves between each pulse. 
For the comb to be useful as a high-precision frequency scale, /cp 
must be stabilised and /ceo should be either stabilised or moni- 
tored. The advantage of the mode-locked laser is that both /ep and 
/ceo are in the radio frequency range so they can be locked (or mon- 
itored) using relatively simple electronics. Such a comb provides a 
series of regular and closely spaced modes, the absolute frequen- 
cies of which are known to an accuracy limited only by that of the 
radio-frequency reference. 
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Figure 1. Laser frequency comb set-up with AAT/UHRF. The comb is fully specified by the repetition frequency, /^p, and carrier envelope offset (CEO) 
frequency, /ceo- The former is measured directly with a photo-diode and counter. The second harmonic generation (SHG) stage frequency-doubles the 
femtosecond-pulsed infrared laser (1.56 \aa) to 780 nm (2f) and also produces some third harmonic light (3f) for the CEO measurement: the high-intensity, 
780 nm light has its fairly narrow bandwidth broadened in a non-linear micro-structured crystal fibre for mixing with the 520 nm (3f) light. The frequency 
difference ('beat note') between them specifies the CEO frequency. Ten out of every 11 comb modes are interferometrically suppressed using a Fabry-Perot 
cavity, increasing the eftective mode separation from 94 MHz to 1 GHz, wide enough to be resolved from each other by the UHRF in its R x 940 k mode. The 
cavity's free spectral range (FSR) is locked to the filtered comb's effective repetition frequency by the servo electronics. The filtered comb light is transported 
to the UHRF pre-slit chamber where it is launched through a diffuser onto the spectrograph slit. 



2.2 Laser frequency comb setup 

The first frequency combs were based on solid-state Ti:Sapphire 
crystals (e.g. |Ma et al.|2004l l which, despite having exceptional fre- 
quency stability, typically require extensive supporting infrastruc- 
ture and careful expert maintenance. In contrast, fibre-based mode- 
locked lasers are low-cost, have 'turn-key' operation and, impor- 
tantly for demonstration at remote astronomical observatories, are 
portable. For these reasons we employ a commercial fibre-laser for 
this experiment using a stabilisation scheme developed at the Uni- 
versity of Western Australia ( ■Locke et al.|2009i ). This laser delivers 
a 100-fs duration pulse at 1560 nm with /^p = 94 MHz repetition 
frequency at an average power of approximately 250 mW. /^p was 
stabilised by detecting its 100 th harmonic on a very fast photodi- 
ode and referencing to a stable microwave signal derived from a 
high quality synthesizer (see Fig.[TJ. 

The laser output is frequency-doubled in a periodically-poled 
lithium niobate (PPLN) crystal, simultaneously producing a red 
comb with 90 mW at 780 nm together with a green comb with 
200 [i-W at the third harmonic, 520 nm. The extraction of the off- 
set frequency is as follows (illustrated in Fig. [TJ. Firstly, the green 
(third harmonic) comb is separated from the red (second harmonic) 
comb and travels a separate path on the optical table. Meanwhile, 
the red comb is frequency-shifted to the green wavelength (520 nm) 
using a unique, highly non-linear optical fibre developed at the Uni- 
versity of Bath ( |Benabid et al.|20()8) . A delay arm in the beam path 
is adjusted so that femtosecond pulses from both combs arrive at 
an avalanche photodiode simultaneously. The difference-frequency 
generated by interference of the green third harmonic comb and 
the frequency-shifted red comb (both at 520 nm) contains the offset 
frequency. 

Although it is possible to feed this /ceo measurement back 
into the laser, thus actively stabilising /ceo, it was not necessary 
for the purposes of our LFC demonstration and IPS measurement. 
This minimized the equipment required at the remote AAT site. Be- 



cause /cEo is not multiplied by n (unlike /cp), it was easily logged 
and monitored with simple electronics during our experiments. In 
principle, this record could be used in post-processing to calibrate 
the comb, though for the purposes of our analysis in Section|4] this 
was not necessary. 

One important characteristic of the frequency comb used here, 
as with all previous tests of LFCs on astronomical spectrographs, 
is that /(.p is too low for the teeth of the comb to be resolved 
from one another by the spectrograph (e.g. [Steinmetz et al.|2008[ 
[WiUcen et al.|2010^ . Indeed, this was one challenge for the imple- 
mentation of 'turn-key' LFCs in astronomy discussed by I Murphy] 
|et al.| ( |2007] >. Our LFCs repetition frequency, /ep = 94 MHz, is 
too low even for UHRF's highest resolution mode with i?=940000, 
i.e. FWHM=0.4GHz. A simple solution employed in previous 
works is to filter out most modes using a Fabry-Perot (FP) cav- 
ity. We chose the FP cavity's free spectral range (FSR) to transmit 
only every 1 1th comb mode, leaving an efi"ective mode separation 
of 1 .03 GHz - wide enough to resolve with UHRF. In designing 
the FP cavity it was firstly coupled to a continuous wave laser for 
end mirror alignment to form a stable cavity. For it to operate with 
the mode-locked laser a further condition must be satisfied, that is, 
every 1 1th (in our case) pulse must map back precisely onto itself, 
requiring a precision in the length of the cavity of order of a few 
wavelengths. One of the mirrors on the FP cavity was set on a fine- 
screw to obtain coarse cavity length tuning, and then the optical 
output was monitored and used to control the cavity length via a 
piezoelectric transducer (PZT) on the mirror. 

The dynamic range of linear response for astronomical CCDs 
is limited to ~40000 photo-electrons, the CCD studied here being 
no exception. Our LFC generally supplied more photons than re- 
quired in short exposures; each comb exposure only needed to be a 
few seconds long. Therefore, individual exposures have very high 
SNR, dominated by the LFC photon statistics, thereby enabling 
constraints on the IPS variations. We are also able to test the stabil- 
ity of the spectrograph on short time-scales using a series of short 
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exposures (see Section [5T| l. By comparison, the usual wavelengthi 
calibration method - a hollow-cathode ThAr emission-line lamp 
- provides only 3 emission lines near 780 nm which provide only 
moderate SNR (~50per pixel) even after several minutes of expo- 
sure. 

Finally, we note that our LFC setup enables a simple method 
for 'flat-fielding' (dividing out the pixel-to-pixel sensitivity varia- 
tions in) the CCD exposures: we simply bypass the FP cavity to 
provide a high-intensity spectrum whose spectral structure is un- 
resolved by UHRF. Again, the usual flat-field calibration method, 
which utilizes a quartz lamp, would have involved averaging many 
exposures lasting several minutes each. 

2.3 UHRF spectrograph setup 

The UHRF is a (grating) cross-dispersed echelle spectrograph 
mounted in the (east) coude room of the AAT (Diego et al. 1995). 
It was selected for our LFC test because its very high resolution 
mode (/?=940000) allows us to use a relatively inexpensive, com- 
mercial LFC with a 'low' (actually typical) /iep. The demands on 
the finesse of the FP cavity mirrors (which was F ~ 300) are also 
reduced because 'only' 10 out of every 1 1 comb modes need to be 
interferometrically suppressed. And while UHRF in this mode only 
allows ~ 1 nm of a single echelle order to be observed in a single 
exposure, a wider bandwidth was not required for measuring pos- 
sible IPS variations. Also, pre-defined echelle and cross-disperser 
angles were already available for observing at 780 nm, so minimal 
time was spent configuring the spectrograph itself. 

The FP-filtered comb light was coupled into a single-mode 
optical fibre for transport to UHRF's pre-slit chamber where it was 
launched in free space through a rotating paper diffuser to over- 
fill the slit in the spectral direction and almost fill it in the spa- 
tial direction. The full UHRF resolving power of i?=940000 was 
achieved by using the smallest slit-width possible, just 40 vim, or 
0.059 arcseconds projected on the sky. The slit had an unvignetted 
length of ~3.4mm, or ~5 arcseconds. The alignment of this sim- 
ple launch setup was checked only by eye using a television cam- 
era system to view the 780 nm light falling on cards placed along 
UHRF's optical axis. The laser light (almost) filled the UHRF col- 
limator symmetrically about its centre. The alignment was not al- 
tered throughout the exposure sequences reported in this paper. 

Once the optical alignment of the launch system was finalised, 
the CCD dewar was manually rotated around the optical axis and 
fixed so that UHRF's spectral direction roughly aligned with the 
CCD rows. This simplified the analysis by ensuring that comb 
spectra from individual CCD rows spanned the entire length of the 
CCD. 



2.4 Charge-coupled device (CCD) details 

The CCD studied here, called 'MITLL3' at the AAT, is chip 2 from 
wafer 67 (i.e. device W76c2) of fabrication run CCID-20 from 
MIT/Lincoln Labs. It has 2048x4096 15xl5-p.m pixels and was 
designed to have a deep depletion layer (i.e. high resistivity), yield- 
ing pixels which are effectively 40 \im thick. 

Several trade-offs were considered in selecting the MITLL3 
chip over another commonly-used chip at the AAT, a back-side il- 
luminated, thinned EEV chip with 2048x4096 13.5xl3.5-vim pix- 
els. The latter is used for precise radial velocity work at the AAT 
(e.g. [Tinney et al.||2003) l which might benefit from knowing the 
IPS variations in that device. Nevertheless, chips from the same 



fabrication run (CCID-20) as the MITLL3 are used at many tele- 
scopes world-wide in many different instruments including, for ex- 
ample, the Ultraviolet and Visual Echelle Spectrograph at the Very 
Large Telescope, Chile ( [Dekker et al.|2000l l. Another consideration 
is that, being thinned and back-side illuminated, the EEV chip's 
IPS variations may be larger and more easily detected by our LFC 
test. Conversely, IPS variations in the deep-depletion MITLL3 chip 
might be harder to detect. 

However, the main factor in selecting the MITLL3 chip was its 
very low fringing level at 780 nm. The much larger fringing level of 
the EEV would have prevented adequate flat-fielding of our comb 
exposures and precluded any measurement of the IPS variations in 
the EEV chip. 

For reference, the MITTL3 CCD pixel size in the spectral di- 
rection of 15 nm translates to 2.63 mA or, at 780 nm, 129 MHz or 
101ms"'. The nominal resolving power of R = 940000 at 780 nm 
translates to a FWHM resolution of 8.3 mA, 409 MHz or 3 19 m s"' . 
The FWHM resolution is sampled by 3.15 pixels. 



2.5 Data reduction 

Before any analysis of IPS variations, only basic data-reduction 
steps were required. The bias level was estimated using the me- 
dian of 51 zero-second CCD exposures and subtracted from each 
CCD exposure. As mentioned in Section [2^ above, flat-field ex- 
posures were obtained by bypassing the FP cavity stage. A master 
flat-field, formed from the median of 3 1 such exposures with high 
photo-electron counts, was used to normalize each LFC exposure 
used in our analysis. Part of a bias-corrected, flat-fielded LFC ex- 
posure is shown in Fig.|2] Note that the spectrum's spatial direction 
is tilted with respect to the CCD columns, a fact which assists our 
subsequent analysis. 

It is important that the error spectrum - the flux uncertainty in 
each CCD pixel - is estimated for our analysis of possible IPS vari- 
ations. The dominant noise source is the Poisson noise of the LFC 
photons; the CCD's read-noise was only ~3 photo-electrons per 
pixel. Thus, the error spectrum varies strongly along the spectral 
direction of any given CCD row. An error image was generated for 
each bias-corrected, flat-fielded LFC exposure by taking the square 
root of the bias-corrected flux image and normalizing by the master 
flat-field image. The error spectrum is shown for a single CCD row 
inFig.|2] 



3 MODELLING THE FREQUENCY COMB CCD 
EXPOSURES 

To understand the extent of IPS variations in the CCD pixels, we 
model the frequency comb flux using non-linear least squares 
minimisation (e.g. |Fisher|1958| l. The reliability of this technique is 
demonstrated using simulated comb CCD recordings where rele- 
vant in the following sections. It is possible that other techniques 
would be effective, perhaps even simpler or faster but, for the 
demonstration purposes of this paper, a full exploration of other 
techniques is not important. The minimisation code is custom- 
written but employs the standard Levenberg-Marquardt method 
with an implementation following the basic structure of that set out 
in|Presset^([T992|. 

We begin by modelling the UHRF instrumental line-shape 
(or instrumental profile) from the 1 -dimensional comb spectrum 
in individual rows of the comb CCD exposures (Section [3.1^ . That 
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Figure 2. Part of a CCD exposure at 780 nm from the AAT/UHRF frequency comb test (top panel) and the electron count spectrum extracted from a single 
CCD row (y = 750; second panel from top, black histogram). The raw CCD image was bias-corrected and flat-fielded before plotting. The spatial direction 
is clearly tilted with respect to the CCD columns. The lower two panels show a zoomed-in portion of the upper two panels. The lowest panel shows a small 
number of the comb modes transmitted through the Fabry-Perot cavity (see Fig.[TJ which have a spacing of ~8pix in the spectral direction (,r coordinate). 
Plotted below the single-row spectrum (black histogram) is the en'or spectrum derived from photon statistics (red line) exaggerated by a factor of 50. 



model is constrained by the assumption that other observed spec- 
trum properties (e.g. average flux, frequency dispersion per CCD 
pixel) vary only slowly across the CCD. This basic model of the 
instrumental line-shape is then fixed (though its parameters remain 
free) and the IPS variations are introduced to the comb model in 
Sectionim 

3.1 Modelling the instrumental line-shape (ILS) 

For our purposes, the instrumental line-shape (ILS) is defined as 
the profile formed in a single CCD row by an infinitely narrow LFC 
mode after passing through the spectrograph. Note that our LFC's 
modes are completely unresolved by UHRF and so each recorded 
comb mode will directly represent the ILS at its observed wave- 
length. The frequency stability of our LFC implies that the spec- 
tral width of each comb mode is ~ 1 MHz, or ~0.8 m s"' at the ob- 
served wavelength of 780 nm. Given the nominal UHRF resolution 
of ~320ms"', the comb modes are very much unresolved. 

For a single CCD row, we initially model the comb flux, F, (in 
electron counts), in CCD pixel ; as the product of an 'envelope in- 
tensity' function, E(x), and a 'comb function', C(x), both of which 
are functions of position in the spectral direction of the CCD, x: 

J^;"+l/2 /^/+l/2 
E{x)C{x)dxx Ei \ C(x)dx. (2) 
1-1/2 Ji-1/2 

The sensitivity of F, to IPS variations is directly linked to the 



fact that the comb function varies quickly on the scale of a single 
CCD pixel, so it is important to integrate the model of the comb 
across the pixel rather than approximate it by taking, e.g. its value 
at the pixel's centre. However, the envelope intensity of the comb, 
E(x) - the variation in the measured peak intensity from mode to 
mode - varies very slowly along the spectrum; the approximation 
in the right-hand-side of equation |2| reflects this. For all the fits 
conducted in this paper, E{x) is simply approximated as a second 
order polynomial. 

Because the LFC modes are much narrower that the UHRF 
resolution, in equation ^ the comb function, C(x), can be written 
as a sum of ILSs from each comb mode, 

C(x) = J^ILS(x,x„). (3) 

The iith mode of the comb transmitted through the Fabry-Perot cav- 
ity occurs at pixel position x„ given by 

X„ = Xo + COr^ + CL)2n^ + . . . , (4) 

with a),, the 'pixel-space repetition frequency'. Equation Q is anal- 
ogous to equation |TJ but, because we model the comb in pixel 
space rather than frequency space, higher-order terms are required, 
characterised by 0)2, 0)3 etc. . In this work we found no need to 
include terms beyond second order in equation l|4| because we 
model such short frequency ranges of comb spectra. Visual inspec- 
tion of Fig.|2]reveals that comb modes fall on approximately every 



© 2012 RAS, MNRAS 000,[T|fTT] 



6 M.T. Murphy et al. 



7^6x10* 




1950 



2000 



2050 



Spectral coordinate, x [pixels] 



Figure 3. Modelling the instrumental line-shape (ILS) from small sections of comb spectra. The plot shows single- and double-Gaussian ILS models (thick 
green and narrow red lines, respectively; the latter lies almost entirely over the former) fitted to part of row y = 750 of a single flat-fielded comb exposure 
(black histogram). The residuals between the two models and the data (again, single-Gaussian in thick green, double-Gaussian in narrow red), normalized to 
the 1-cr errors, are plotted below the spectrum relative to the ±3-o" levels (dotted blue lines). Although the two fits are almost indistinguishable by eye, the 
residual spectra show that the double-Gaussian ILS model is strongly preferred. 



8th pixel - i.e. oj^ ~ 8pix - without any obvious non-linearity, 
i.e. na)2<^ <ji>r- Indeed, in our detailed modelling below, coi is typi- 
cally ~ -1 X 10"* with -50 :S n < 50. 

Visual inspection of the example comb spectrum recorded in 
a single CCD row (see Fig. |2j reveals UHRF's ILS to be slightly 
asymmetric. Our initial guess for the ILS was therefore a double- 
Gaussian model: a central Gaussian flanked by a weaker second 
one with the same width. That is, 

M 

ILS(x, x„) = ^ A„, exp [-{x - A„, - x„f jla]^^ (5) 

m=\ 

with M = 2, Al =0 and the ratio of Gaussian amplitudes, AijAi, 
normalized so that J^^ ILS(x, x„)dx = 1. 

Figure [3] shows the model specified in equations ([2]|5]( fitted 
to a small section of a single row of one comb exposure. Also 
plotted is the same model but with a single Gaussian ILS instead 
of a double-Gaussian (i.e. M = i instead of 2). The two models 
are not readily distinguishable by eye. However, the residual spec- 
tra clearly show that the double-Gaussian model provides a closer 
match to the data. This is also reflected in the values of psr de- 
gree of freedom, v, for the two fits: xi. = 29.5 and 2.7 for the single 
and double-Gaussian models respectively. We discuss why the lat- 
ter value is noticeably greater than unity in Section [4~2l The resid- 
ual spectrum for the double-Gaussian model, also plotted in Fig.|3] 
does not appear to show any features which repeat on the scale of 
the comb mode structure, unlike that for the single-Gaussian model. 
We attempted to include a third Gaussian (i.e. M = 3) in the ILS but 
its best-fitted amplitude was typically very small (A3/A1 < 0.02) 
and, upon fitting different CCD rows and using different CCD ex- 
posures, its separation from the main Gaussian component was 
sporadic and very weakly constrained. This is consistent with the 
lack of periodic structure in the residuals from the double-Gaussian 
model in Fig. [3] Also, including the third Gaussian barely reduced 
xl in most cases; when fitting some CCD rows from some expo- 
sures it actually increased. That is, there is no clear evidence for 
additional structure in the ILS beyond the double-Gaussian model. 

[Diego et al"] {T995| measured the UHRF ILS using a He-Ne 
laser. From their figure 10, their measured ILS can be described 
as comprising two main Gaussian-like components separated by 
-0.25 km s"' with similar widths of FWHM ^ 0.28 km s"' and 
with one being a; 20 per cent the height of the other. This accords 
well with the ILS model derived from the comb spectra: in the 
simple fits above, the Gaussian components were separated by 
a;0.23kms"', they had FWHM a 0.30km s"' and the weaker one 
was 20 per cent the height of the stronger one. Thus, the double- 



Gaussian model of the UHRF ILS seems well justified, both from 
the comb spectra themselves and the He-Ne spectra of |Diego et al.| 
( fT995) . 

When fitting different parts of the same CCD exposure stud- 
ied in Figs.[2]&[3] and the others taken over the course of our LFC 
test, we found no significant variations in the ILS structure. That 
is, no further Gaussians seem to be required. Nor do we observe 
large variations in the parameters of the two Gaussians constituting 
our preferred ILS model; when modelling a 2-dimensional coinb 
spectrum in Section |4] we explicitly test for this possibility. How- 
ever, throughout this rest of this work we adopt the basic double- 
Gaussian model of the ILS with free parameters as specified in 
equation l|5j. 

Finally, a by-product of the above fits is the statistical uncer- 
tainty on xo in equation the zero-point of the comb function. 
For example, the double-Gaussian fit to the 150 spectral pixels in 
Fig.jsjretums a 1-cr uncertainty on xo of 9.9 x 10"^ pix, or «1 ms"'. 
In the absence of intra-pixel sensitivity variations (which we seek 
to measure in this work), this would correspond to the precision 
with which such short sections of the comb spectrum could be 
wavelength-calibrated. 



3.2 Modelling CCD intra-pixel sensitivity (IPS) variations 

In our comb test with UHRF we sought to detect or limit any IPS 
variations only in the spectral direction across the CCD pixels. 
This is possible because the flux measured in a given pixel will 
depend on the IPS variations most strongly when the comb flux 
changes substantially across a pixel. From the example comb expo- 
sure shown in Fig.|2]it is clear that the flux variations in the spatial 
direction are much more slowly varying compared with those in the 
spectral direction. Therefore, we do not attempt to model or detect 
IPS variations in the spatial direction across pixels. 

Thus, when fitting the comb flux in a single CCD row, or part 
thereof, we insert a term for the IPS map, IPS,(x), in equation |2]l 
to obtain 

w+1/2 

F, « Ei C(x)IPS(x)dx, (6) 

Ji-l/2 

where, as discussed above, IPS(x) is the same for all or IPS,(x) = 
IPS(x). 

Our aim here is to measure the average IPS variations, not the 
IPS map for individual pixels. The latter possibility is discussed 
in Section|5] We therefore make two important assumptions in the 
modelling that follows: 
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(i) The IPS variations in all CCD pixels under consideration are 
assumed to be the same. While previous laboratory studies have 
shown that IPS variations are different in different pixels, they also 
demonstrate a high degree of correlation across large portions of 
the CCD (e.g. |Toyozumi & Ashley|2005) . 

(ii) The IPS variations are modelled as a symmetric function 
about the centre of the pixel (in the spectral direction). We report 
results here from a Gaussian function with variable width, trun- 
cated at the pixel edges and a base fixed at zero with no slope. For 
Gaussian widths larger than FWHM ~ 0.5 pixels, which we find 
are always preferred for the UHRF CCD, a variable base would 
be closely degenerate with the width parameter. Similar functions 
- e.g. symmetric linear slopes away from a central IPS peak (or 
trough) to the pixel edges - give very similar results. It is also im- 
portant to realise that a variable baseline slope would be closely 
degenerate with the parameters in equation |4]( which characterize 
the position of the comb modes. That is, a general slope in the IPS 
map of the fitted pixels is not detectable in our work. We discuss 
this further in Section l?!2l 

We will see in Section [42] that it is only by fitting many hun- 
dreds of pixels (all assumed to have the same IPS map) that the 
degeneracies in our fitted model are reduced to a level where the 
Gaussian width of the IPS model can be constrained. For exam- 
ple, it is the small asymmetry in the UHRF ILS which allows us 
to constrain a symmetric IPS map, as we have chosen to use here. 
Nevertheless, this proves sufficient for measuring IPS variations in 
the MITLL3 chip using the assumptions above. As we demonstrate 
in Section[5] it may be possible to measure non-symmetric average 
IPS maps in future using additional information. 

To summarize the above discussion, we may write our model 
of the IPS map for all pixels as a function of x (where we note that 
X = is the middle of a pixel) as 

IPS,(x) = IPS(x) = B exp (-x-/2o-?ps) , (7) 

where the amplitude, i?, is normalized such that J IPS(x)rfx = I. 
Here, o"ips, measured in units of pixels, is the IPS Gaussian width 
parameter which we measure. For example, crips » 1 pix indicates 
a flat average IPS map whereas crips < 1 pix indicates a strongly 
variable average IPS map. 

Finally, note that the form of equation ^ is convenient be- 
cause the integrand is the sum of a product of Gaussian functions, 
thereby allowing trivial numerical integration using the well-known 
Gamma function (e.g. jPress et al.|1992^ . 



4 AVERAGE INTRA-PIXEL SENSITIVITY MAP FROM A 
SINGLE COMB EXPOSURE 

To measure the IPS variations using a single CCD exposure we as- 
sume that they are the same in all pixels and utilise the fact that each 
pixel (i.e. "copy" of the IPS map) is probed with a different flux dis- 
tribution. By using a well-constrained model of the flux distribution 
across many pixels, crips - the Gaussian width of the average intra- 
pixel sensitivity (IPS) map in our model - can be measured. From 
preliminary fits to ~200-pixel sections of single CCD rows, it be- 
came clear that cTips was covariant (i.e. somewhat degenerate) with 
several other parameters in our comb and instrumental line-shape 
models, e.g. ctils, the width of the two Gaussian ILS components. 
Rather than attempt to fit many more pixels in a single CCD row, 
which may necessitate higher-order representations of the intensity 



envelope function, £(x), and comb function, C(x), we instead uti- 
lized the full 2-dimensional (2D) comb information in our UHRF 
CCD exposures. 



4.1 2D modelling of a single comb exposure 

The example CCD exposure in Fig.|2]illustrates how the comb sig- 
nal occupies more than 100 pixels in the spatial direction. That is, 
there are many "copies" of almost the same comb signal available 
for constraining crips. This means that, under our assumption that 
all CCD pixels have the same IPS variations, we can fit many CCD 
rows simultaneously by adding a very small number of additional 
free parameters to the model in equation (j6j. 

Firstly, each of the terms of the comb function in equation |4]( 
now becomes a function of CCD row, j. That is, the nth mode in 
row j can be written as 

x{ = xo(l + aij + a2f) + Wr«(l + b\j) + (jj2n^ (8) 

where we have now explicitly truncated the series at second order 
in n and limited the order in j in each term. Figure[2]clearly shows 
that the spatial direction of the spectrograph is tilted with respect 
to the CCD columns, i.e. ai > to account for this tilt. But after 
exploring many 2D fits to various comb exposures and portions of 
the CCD, we found no evidence for distortions in this tilt, or for 
variations in the pixel-space repetition frequency with row. That 
is, both ai and b\ were consistent with zero. Therefore, for all fits 
reported here we fixed ai and bi to zero. 

Similarly, each of the parameters of the instrumental line- 
shape (ILS) could vary with j. By appropriately altering equation 
ijsj to include polynomial dependence on j for all free parameters 
(e.g. (Tils), w6 explored various fits to search for evidence of this. 
However, we found no significant variations in the ILS in the spatial 
direction. 

Finally, note that while the intensity envelope function, E(x), 
clearly varies (slowly) with CCD row, y, it will not be strongly co- 
variant with any IPS variations, the measurement of which is our 
main aim. Therefore, we simply allowed the parameters of E{x) to 
be independently determined for each j. 



4.2 Main results: limits on IPS variations from individual 
exposures 

We searched for IPS variations in our 10 highest SNR (20-s 
long) exposures. In each exposure we concentrated on fitting a 
50x320 pix (rows x columns) region near the centre (highest SNR) 
of the comb signal. This represents about half the comb's spatial 
extent and about 10 per cent of its spectral extent on the CCD. Pre- 
liminary fits showed that some portions of the fitted area were not 
well flat-fielded. From inspection of the raw images and the mas- 
ter flat-field frame, we determined that dust particles on the UHRF 
optics were most likely resulting in small flat-field residuals which 
varied slightly from exposure to exposure or at least between the 
comb and flat-field exposures. The portions of the fitted area were 
the same in the 10 exposures we studied in detail. Therefore, we 
simply masked these portions out of om minimization analysis. 
We note in passing that 2D fits to different regions of the CCD did 
not reveal any obvious periodic CCD defects like those identified 
by |Wilken et"aLl ( |20T0) . 

Figure|4]shows part of a 2D fit to a representative exposure. In 
general, the fit is clearly a reasonable representation of the recorded 
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Figure 4. Part of a 2D fit to a single comb exposure. The flat-fielded spectra from every 6th fitted row are shown as black histograms overlaid by the model 
(green line) which was fitted to 50 rows and 320 pixels in the spectral direction. The normalized residuals are plotted below the spectrum of each row relative 
to the ±3-o" levels (dotted blue lines). Note that some sections of residuals are all zero; these parts of the CCD have been masked due to problems with 
flat-fielding (see text). Note that the tilt of the spatial direction across the CCD columns is evident. 



data. The;^'^ per degree of freedom, v, was similar for all 10 expo- 
sures, xi ~ 2.9 which, while not «1, is not unreasonable, espe- 
cially in the possible presence of additional, low-level flat-fielding 
errors. The residual spectra generally do not show coherent struc- 
tures which repeat over ~8 pixel ranges, i.e. the transmitted comb 
mode spacing. Figure|4]clearly shows some exceptions to this gen- 
eral statement. However, we did not observe any obvious correla- 
tion between the positions of these exceptions in different expo- 
sures. They may very well result from IPS variations in small clus- 
ters of pixels, but we cannot quantitatively address that possibility 
with the data in hand. Again, we attempt here to constrain a model 
of the average IPS variations for the whole fitted area. 

For all of the 10 individual exposures, simply including iTips 
as a free parameter in the fit yielded best-fit values cTips > 2 pix in 
all cases. This was independent of the starting guess for ctips (and 
other parameters in the model), which was generally ~0.5 pix. In 
most cases crips grew from the starting guess value to 100 pix, the 
arbitrary limit imposed in our code. And in all cases, the best fit 
value was entirely consistent with crips = oo. That is, under the 
assumptions of our IPS model, we do not detect significant average 
IPS variations. 

To ensure that our;^'^ minimization code was functioning cor- 
rectly we fitted simulated versions of the 10 comb spectra. The real 
comb flux was replaced with the best-fit model (without IPS varia- 
tions) with Gaussian noise added which was derived using the real 
flux error estimates. The synthetic realisations were fitted in the 
same way, and with the same free parameters, as the real spectra. 
Figure [5] illustrates one aspect of the comparison between the re- 
sults from one comb exposure and the average results from 25 sim- 
ulated versions of it. By holding crips fixed at certain values and 
plotting the deviation in;^'^ from its minimum, Ax^, we gain a more 



accurate understanding of the minimum value of crips allowed for 
each exposure. Because xl is not close to unity for the real expo- 
sures, but it is (by construction) for the simulations, we can test the 
veracity of the usual assumption that confidence intervals for trips 
may be derived after normalizing by CtJ,^i„)''". the square root 
of the minimum Under this normalization. Fig. [s] shows very 
similar results for the real and simulated spectra. 

For the particular exposure studied in Fig.|5]we can rule out 
(Tips < 1-0 pix (< 0.5 pix) at the 68 per cent (95 per cent) signif- 
icance level. These constraints on crjps were uniformly replicated 
for the same fitted region across the 10 exposures we studied in de- 
tail. By co-adding the plots of Ax~ /(xl^i^y^^ versus crips for all 10 
exposures, and recalling that we fitted the same region of the CCD 
in all cases, we form a joint constraint on the Gaussian width of our 
assumed IPS variation map: crips > 1.2 pix at 99.7 per cent (3-cr) 
confidence. This value of a"ips corresponds to variations in the IPS 
across the spectral direction of the pixel of < 8 per cent. 

Instead of a decrease in the IPS from the pixel centre to its 
edges, we also explored an IPS map with a minimum at the pixel 
centre by adding unity to equation and allowing S < 1. We 
find very similar results to those above: trips is constrained to be 
> 1.3 pix at 99.7 per cent confidence. That is, we rule out IPS vari- 
ations of > 7 per cent at 3-tr significance under the assumptions of 
our modified (up-side down Gaussian) model. 

Finally, to demonstrate that detecting non-zero IPS variations 
is possible using our^^ minimization technique (under the assump- 
tions of our 2D modelling), we again fit simulated versions of the 
real comb spectra, this time with finite values of trips. The SNR of 
the simulated spectra again matched the corresponding real spectra. 
In all simulations the best-fit value of tr ips broadly matched the in- 
put value. This was simple to assess for input values trips ^1-5 pix 
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Figure 5. Results from 2D fit to the central part (50x320 pixels) of a single 
flat-fielded comb exposure. The plot represents the increase in from its 
minimum value, . (i.e. Ay^ = - . ) as a function of the Gaussian 

^min ^ /V /V ^min' 

width of the intra-pixel sensitivity (IPS) variations, crips. is normalized 
by -yjxyjni^ "^^^^^^vmin minimum value ofx^ per degree of freedom, 
xl, so that confidence intervals can be calculated directly from the plot, crips 
is unbounded above in the fit to the real comb exposure (blue curve with cir- 
cular points). That is, we cannot rule out the null hypothesis of zero average 
IPS variations in that exposure. However, trips 0.5 pix is ruled out at 
95 per cent confidence {Ax^/ ^vfmin ~ A fit to a simulated version of 
the same comb exposure with a flat IPS map inserted yield very similar re- 
sults (red curve with triangular points). The green curve with square points 
shows the results from a fit to the same simulation in which a non-flat IPS 
map was used, i.e. crips = 1 pix. The definite minimum demonstrates that 
non-flat IPS variations could have been detected given the assumptions in 
our model of the comb exposures. 



but judging formal statistical consistency between the input and 
best-fit values was more difficult when ctips was larger. The ultimate 
reason for this is that the IPS map is very flat for large (Tips. This is 
illustrated in Fig. [5] which shows the variation in;^'^ with fixed fitted 
values of ctips for a single realisation with an input ctips = 1 . In this 
example there is a clear minimum in^- near the input value of crips . 
However, it is also clear that^-'s sensitivity to iTips diminishes with 
increasing crips, i.e. the^- curve is not parabolic at iTips > 1. Thus, 
the formal statistical error on cTips derived from the covariance ma- 
trix in our minimization scheme will not be reliable for large 
(Tips. This problem clearly also applies to the constraints on crips 
derived from our real comb exposures above. This is the reason we 
used the;f- curve for the fit to each exposure, rather than the formal 
error from the covariance matrix, to derive confidence limits on any 
IPS variations. 

To summarize the results above, we have demonstrated that 
any symmetric component of the average intra-pixel sensitivity 
(IPS) variations can be detected in a single LFC exposure provided 
that many hundreds of pixels can be modelled simultaneously. The 
AAT/MITLL3 chip's average IPS map deviates from being flat by 
< 8 per cent in a symmetric manner about the pixel centre in the di- 
rection along CCD rows. However, under our model assumptions, 
we cannot measure or rule out a significant asymmetric component 
to the average IPS map. For example, the IPS map may slope from 
one side of the pixel to the other, but we cannot detect this by mod- 
elling a single exposure. The simple reason for this is that such a 
slope is closely degenerate with the measured position of the comb 
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Figure 6. Shifts between comb spectra on short time-scales. Plotted are 
small sections of a single row of three 1-s long flat-fielded comb spectra 
(traced by spline curves to between distinguish them). The exposures were 
taken approximately 1 min apail and all three had the same repetition fre- 
quency (/rep) and carrier envelope offset frequency (/ceo)> yet they show 
clear, stochastic shifts of ~0.4 pixels. 



with respect to the CCD pixel grid. We return to this point in Sec- 
tion l5.2l In the next Section we discuss how a series of olfset comb 
exposures might allow that degeneracy to be removed. 



5 INTRA-PIXEL SENSITIVITY MAP FROM A SERIES 
OF OFFSET COMB EXPOSURES 

5.1 Observed series 

With a view to detecting any very small average IPS variations, 
the 10x20-s comb exposures with high SNR utilized in the pre- 
vious section were taken in series with small offsets in the LFC 
signal applied between exposures. The offsets were generated by 
slightly altering the repetition frequency of the mode-locked laser, 
,/rep, which was observed to be stable between exposures. The off- 
set in /rep between each exposure was ~15Hz so that the comb 
modes at 780 nm would shift by ~50ms"', or 0.5 spectral pixels. 
The carrier-envelope offset frequency, /ceo, was also monitored 
and fluctuations were negligible (< 10 ms"') during the series of 
exposures. 

However, upon inspecting the series of offset comb exposures, 
it was apparent that the observed shifts in the comb modes did not 
accurately reflect those expected from the offsets applied to /ep 
between successive exposures. To investigate this we took another 
series of 10 short (1-s) exposures with no offsets in /rep applied 
between them. That is, including the CCD read-out time of each 
exposure, a 1-s exposure was taken approximately every minute 
and, since neither /ep nor /ceo changed appreciably during the 
10 min series, we expected no shifts between the comb modes. 
However, we discovered fairly large stochastic shifts between suc- 
cessive comb exposures. These were typically <0.2 pixels but some 
were as large as ~0.4 pixels. An example from 3 successive 1-s ex- 
posures is shown in Fig. [6] Furthermore, no two successive expo- 
sures had zero shift between them, indicating that the time-scale for 
these shifts is substantially less than 1 min. 

The origin of these shifts is not clear from our comb spectra, 
but the fast time-scale suggests that a purely thermal drift within 
UHRF is an unlikely culprit and that a more likely explanation is 
small mechanical oscillations or jumps/relaxations in the spectro- 
graph. UHRF's main components and camera are mounted on a 
large (~3-m), vibration-damped and isolated bench. The slit assem- 
bly is mounted on the same bench. However, the pre-slit area upon 
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which our LFC launch setup was mounted (see Fig. [TJ is, at least 
in principle, decoupled from the spectrograph bench. Thus, oscil- 
lations induced by moving components in the pre-slit area (e.g. the 
rotating difFuser) should not cause the slit itself to oscillate. How- 
ever, it remains possible the other equipment in the pre-slit area 
might transmit some vibrations to the slit assembly, which may 
have been responsible for the short time-scale, stochastic shifts we 
observe. If such oscillations/vibrations did not originate from our 
equipment, they may also affect astronomical observations with 
UHRF; however, they are unlikely to have been detectable in UHRF 
before because of the much longer ThAr exposure times required 
(i.e. several minutes). Observations with the University College 
London Echelle Spectrograph, UCLES, which uses the same coude 
room and slit assembly as UHRF, might also be affected. We would 
be interested to learn whether users of UCLES have identified such 
problems before, but to our knowledge they have not been reported 
in the literature. 

Unfortunately, these stochastic shifts between exposures mean 
that we cannot reliably model the observed series of offset comb ex- 
posures to constrain IPS variations. One should consider whether 
the analysis of the individual 20-s exposures in Section |4] was af- 
fected by this instability. The effect of quickly varying shifts in 
the comb position, averaged over a comparably long exposure, 
should be to increase the observed comb mode width, or ctils as 
parametrized in our comb models. We do not observe a larger ctils 
for the 10x20-s comb exposures compared to the lOxl-s ones: 
(ciLs) = 1-632 compared to (ctils) = 1.629 pix, respectively. 
However, the root-mean-square (RMS) variations in these quan- 
tities - 0.006 and 0.01 1 pix, respectively - may be too large to de- 
tect the increase in width expected: if random comb shifts have 
RMS ~ 0.2 pix then ctils would broaden by just ~0.0I pix. Even if 
this effect were present in the 20-s comb exposures, our model of 
the comb flux across the CCD pixels is still likely to be adequate 
and so the constraints on IPS variations we derive are unlikely to 
be strongly affected. 



5.2 Simulated series 

While the short time-scale instabilities in UHRF prevent us from 
constraining the average IPS variations using a real series of offset 
comb exposures, we can nevertheless demonstrate this possibility 
using simulated comb spectra. We simulated 10 spectra, each with 
a single CCD row 320 pixels long. The comb signal was generated 
using the best-fitting parameters derived from row 750 of the CCD 
exposure studied in Section |4.2| (see Figs. [2] & [3j, including the 
double-Gaussian instrumental line-shape. Gaussian random noise 
was added to each spectrum independently, with the average SNR 
set to the same value as in the real CCD exposure. However, we in- 
troduced 0.8-pixel offsets between the comb functions for succes- 
sive spectra. We then simultaneously fitted comb functions to all 10 
spectra, the key assumption being that the relative offsets between 
spectra were known; if the LFC system is absolutely calibrated then 
this will be true to very high precision. 

Figure |7] shows the results of the fits to the simulated spectra. 
Similar to the simulations in the previous section, the input IPS 
model was a Gaussian with ctips = 0.4 pix centred at x = 0. That is, 
the Gaussian was centred at the middle of the pixel and cTips was a 
free parameter. Note that the ILS and envelope function parameters 
were also free in this fit. Figure [7] shows that this symmetric IPS 
model is recovered with a relatively small uncertainty in crips of 
0.03 pix. This demonstrates that a series of offset comb spectra can 
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Figure 7. Recovery of average intra-pixel sensitivity (IPS) map from se- 
ries of offset comb spectra. Ten single-row, 320-pixel comb spectra were 
simulated with the same parameters (similar to those of our real comb ex- 
posures), but with offsets of 0.8 pixels between their respective comb func- 
tions and different noise realisations. The blue (symmetric about jr = 0) 
shaded region shows the l-o" uncertainty range in the measured IPS map 
(solid curve) when only the Gaussian width of the IPS model is a fitted 
free parameter The red (asymmetric) region and curves represent the re- 
sults when both the Gaussian IPS width and its position along x are free 
parameters and where it is assumed that the absolute position of the comb 
with respect to the CCD pixel grid is known. In both cases the input IPS 
model is shown as a black dotted curve. 



be used to characterise the symmetric component of the average 
IPS map. 

The main reason a series of offset spectra can strongly con- 
strain the IPS variations, even with fairly few pixels (e.g. a single 
CCD row) assumed to have the same IPS map, is that the comb flux 
distribution across each pixel is changed by a known amount when 
the comb is offset by a known amount. However, as with the 2D 
analysis of a single exposure in Section|4] any asymmetric compo- 
nent of the IPS map is highly degenerate with the overall position 
of the comb. That is, without knowing where the comb modes fall 
with respect to the CCD pixel grid, asymmetries in the average IPS 
map cannot be reliably determined. If some additional information 
could be gained about where the comb light is striking pixels, per- 
haps by adjusting a physical mask in front of the CCD, then the 
same fitting technique could be used to recover asymmetries in the 
IPS map. 

As an example of this, we simulated another ten 320-pixel, 
single-row comb spectra in the same way as above but with an 
asymmetric input IPS map, i.e. a Gaussian, as before, but shifted 
with respect to the pixel centre. It was then assumed that the ab- 
solute position of the comb modes with respect to the pixel grid 
was known. Figure |7] shows that the fitting technique indeed recov- 
ers the asymmetry in the input IPS model with high accuracy. Much 
higher precision measurements of the average IPS variations should 
be possible if the instrumental line-shape (ILS) is also known ac- 
curately a priori (or measured independently) and not simultane- 
ously constrained by the comb spectra. With a large enough number 
of high-SNR exposures with very small, known offsets, and with 
known comb-mode-CCD-grid positioning and ILS information, it 
should be possible to measure the IPS map of every illuminated 
CCD pixel. 
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6 SUMMARY 

We reported the successful test of a laser frequency comb (LFC) 
system on the Ultra High Resolution Facility (UHRF) spectro- 
graph at the Anglo-Australian Telescope (AAT). Figure [T] shows 
the experimental setup used. The high spectral resolution of UHRF 
{R X 940000) allowed a commercial fibre laser with a fairly typical 
repetition frequency (/rep) to be used. 

A Fabry-Perot cavity was used to widen the effective sepa- 
ration between modes to ~lGHz (by suppressing 10 of every 11 
modes) so they could be resolved apart by UHRF. The LFC spectra 
were recorded with the 'MITLL3' CCD chip at s;780nm in a sin- 
gle echelle order of UHRF covering a 10 A. An example exposure 
is shown in Fig.|2] 

Our main goal in the LFC test was to demonstrate that LFC 
spectra could be used to characterize the average intra-pixel sen- 
sitivity (IPS) map of astronomical CCDs in situ. This meant that 
actively self-calibrating the LFC (i.e. controlling /ep and /ceo, not 
just measuring/monitoring them) was not necessary, which mini- 
mized the equipment required on-site. 

We demonstrated the measurement of IPS variations using a 
single comb exposure by modelling the comb flux in 50x320 pixels 
in the spatial and spectral directions, respectively - see Figs.|4]&|5] 
For a simple Gaussian model of the average IPS variations where 
the sensitivity falls off symmetrically about the pixel centre, the IPS 
map of the AAT/MITLL3 chip deviates from flat by < 8 per cent. 
A similar constraint on an inverted map (i.e. sensitivity increasing 
towards the pixel edges) was also obtained. 

The instrumental line-shape (ILS) was also modelled in this 
process, and we found that a double-Gaussian ILS was sufficient to 
fit our comb spectra - see Fig. [3] This ILS is very similar to that 
found when commissioning UHRF ( [Diego et al.|1995| >. 

By slightly adjusting the repetition frequency between comb 
exposures (by changing the length of the fibre laser's free-space 
section), we obtained a series of exposures with known offsets be- 
tween their respective comb signals along the spectral direction of 
the CCD. In this way we intended to more strongly constrain the 
IPS variations in a single CCD row. However, a series of shorter 
exposures without offsets between them revealed short time-scale 
(<1 min) random shifts between exposures of up to «0.4 pixels (or 
a40ms"') - see Fig. [6] This precludes a thorough analysis of the 
series of offset comb exposures for constraining IPS variations. 
Nevertheless, we used simulated spectra whose parameters closely 
matched our real spectra to demonstrate that, with a stable spec- 
trograph, such a series of comb exposures should allow symmetric 
departures from a flat average IPS map of a single CCD row (320 
pixels) to be strongly constrained - see Fig.|7] 

Despite demonstrating that the symmetric component of av- 
erage IPS variations can be constrained using LFCs, it is clear 
that, without knowing a priori where the comb modes fall with 
respect to the CCD grid, any asymmetric component of the IPS 
map cannot be reliably measured. This could be an important as- 
pect of IPS variations to measure for some applications, especially 
those where long-term, effectively absolute calibration is required 
(e.g. the real-time observation of the drift in Lyman-ff forest lines 
in quasar spectra). In future, perhaps the relationship between the 
physical position of the comb modes and the CCD grid might be 
established through some kind of physical CCD masking proce- 
dure, or by 'switching off' some CCD pixels and obtaining a series 
of offset comb exposures. Exploring such possibilities may be re- 
quired to fully realise absolute wavelength calibration using CCDs 
in astronomy. 
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